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SUMMARY 

There  are  several  bi-directional  situations  where  it  is 
required  to  obtain  a measure  of  correlation.  The  question  has 
been  raised  often  that  the  bi-directional  correlation  coefficients 
so  far  known  are  not  scale  invariant  even  asymptotically. 

Following  Cox’s  procedure  (Cox,  1975),  we  introduce  a new 
correlation  coefficient  which  has  the  desirable  property  of 
being  invariant  under  scale  for  large  samples  with  von  Mises 
marginals.  We  obtain  its  asymptotic  distribution  under  the 
hypothesis  of  independence.  We  examine  its  properties,  and 
give  a numerical  example. 
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A.  D.  BLOSE 
Technical  Information 


1 . Cox's  Procedure 


Cox  (1975)  derived  tests  for  "von  Miseness"  for 
circular  observations  using  the  exponential  distribution 
with  probability  density  function  proportional  to 

exp(a^  cos  x +0^  sin  x + a.^  cos  2x  + 02  sin  2x)  . 

If  X , . . . ,Xn  is  a random  sample,  then  the  agreement  with 

the  von  Mises  density  = B 2 = ® is  thus  tested  from  the 

conditional  distribution  of  (7  cos  2X_.  , T.  sin  2X  J given 

(7,  cos  X . , L sin  X . ) . 

D D 

Let  (X, ,Y, (X  , Y ) be  n independent  observations 
11  n n 

on  a torus.  Mardia  (1975)  introduced  the  correlation  coefficient 
of  (Xi,Yi)  , i = 1, . . . , n as 

2 2 

max(D+,D_)/{  U -Rj_)  (1  -R2M  , (1.1) 

where 

n2D  = {Y  cos(X*  ± Y*)  - nR.R_}2  + f Y sin(X*iY*)}2  . 

*1*1  11  12  ix 

X*  = (X^  -x0)  mod  2tt  and  Y*  = (Y^  -YQ)  mod  2n  , and 

where  R^  and  i*2  are  the  mean  resultant  lengths  and 

X_,  Y-  are  the  mean  directions  of  X.  and  Y.  respectively. 

0 0 11 

(Downs  and Eif ler  (1975)  show  the  equivalence  of  this  correlation 
with  a correlation  coefficient  introduced  by  Downs  (1974)  in 


3 


i 


\ nt.ya  & ( 


a different  set  up) . Downs  and  Eifler  (1975)  point  out  that 
the  distribution  of  (1.1)  will  in  general  depend  on  the 
concentration  parameter  of  the  marginals.  (For  a summary  of 
various  circular  correlation  coefficients,  see  Puri  and  Rao 
(1976))  . 

It  is  shown  that  Cox's  procedure  (Cox,  1975)  provides  a new 
correlation  coefficient  which  is  Scale  invariant  for  large 
samples  with  von  Mises  marginals. 

Consider  the  probability  density  function  introduced  by 
Mardia  (1975)  as 

C exp[x^  cos  (x-y)  + *2  cos  (y-v)  + a cos  x cos  y + b sin  x cos  y 
+ c cos  x sin  y + d sin  x siny]  , 0$x,  y s 2rr  . 

Following  Cox's  procedure,  to  test  the  hypothesis  of 
independence  (i.e.  a=b  = c = d = 0)  with  von  Mises  marginals, 
we  should  consider  the  distribution  of 

nU'=  (E  cos  X.  cos  Y.  , E cos  X.  sinY,  , r.  sinX,  cos  Y.,  I sinX.  sin  Y. 

~ l l l l l l ii 

given 

nT'=  (ycoaX.,  7 sinX^,  r cos  Y.,  EsinY^ 

To  obtain  a correlation  coefficient,  a suitable  function 
of  U(T  is  considered  in  section  2.  Its  asymptotic 


distribution  under  the  null  hypothesis  is  considered  in 
section  3.  The  properties  of  the  resulting  coefficient  is 
examined  in  section  4 together  with  its  spherical  extension- 
In  section  5,  a numerical  example  is  given. 

2.  An  Invariant  Function 

Consider  the  hypothesis  H of  independence:  a=b  = c = d=  0 . 

Under  H , let 

Cov(u,T)  = ~ E = ~ (li11  ^12)  - 

—21  ~22' 

Set 

— 1 • 2 = -11  " -12-22  -21 

and  denote 

,(X.Y)  S^tC-E^-Sj^fT-E^)],  (2.11 

where  E denotes  the  expectation  under  H , and 

X'  = (X. X ) and  Y ' = (Y, . . . . ,Y  ) . 

— In  —In 

A A • A • 

Let  now  E#  y (U)  and  E (T)  denote  the  maximum 

-w  h — H ~~ 

likelihood  estimators  of  E,  E„(U)  and  E„(T)  respectively. 

It  may  be  noted  th?.t  these  estimators  are  the  functions  of 
the  maximum  likelihood  estimators  of  |i,  v,  h , and  *2  where 

A A A 1 A 1 

u = . v = YQ  * x ^ = A (R^)  , and  x 2 = A (2.2) 

in  the  notation  of  Mardia  (1972,  Chapter  5) . Then  the 
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resulting  value  of  g(X,Y)  is 


i 


XY 


CH-‘eh(2»  -ii2 


(2.3) 


The  rationale  for  (2.1)  is  obvious  considering  the  conditional 
distribution  of  U|T  when  U and  T are  jointly  normal.  Note 
that  (2.2)  and  (2.3)  are  invariant  under  separate  linear 
non-singular  transformations  on  U and  T . 

We  now  show  that  this  procedure  when  applied  to  the 
bivariate  case  with 

nU  = ? X.Y.  , nT'=  (TX.,  TX2.  TY . , TY2) 

~ . , i i — i i i i 

i=l 

2 

leads  to  r where  r is  the  ordinary  correlation  coefficient. 

We  assume  that  under  H , X and  Y are  independent 
2 2 

N(u»ct^)  and  Nlu.o^)  respectively.  Then  it  can  be  easily 

2 2 -1 

shown  that  £1>2  = o jO 2 , I12~22  = and  (2.1) 

reduces  to 

(U  - UV  - v (X  - u)  -u(Y  - v)  )2/o2c2  i (2.6) 

n n 

where  X = T X./n  and  Y = T.  Y./n  . Using  the  maximum 

i=1  1 i=1  2 2 
likelihood  estimators  for  p,  v,  and  ct2  under  H , 

2 

(2.4)  reduces  to  r . Hence  Cox's  procedure  does  lead  to  a 
suitable  correlation  for  the  linear  case. 
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3.  The  Correlation  Coefficient 
Without  any  loss  of  generality,  let 


X = Y = ° 
o o 


(3.1) 


Consequently,  yi  = v = 0 . 
It  can  be  shown  that 


7l  us*  0V* 


/ 


A.  1 A A*  ^2^*  A * * 

^ diag{  (l  + a2)(l  + a2)  - 4a  a , (1  + a2)  (1  - a2)  . 


(1 -o2)  (1 -a*),  (1  -o2>  (1  +a*)]  , (3.2) 


~ 22 


diagfa^  a2>  a*,  a*)  , and 


—12 


aa* 


1' 


0, 

0, 

0, 


aa* 

0 

0 


(3.3) 


where 


2al  - 

A A 2 A 

1 + a2  - 2a  , 2a 2 = 1 - a2 

(3.4) 

2a*  « 

* ★ * *2  ^ ★ 

l+a2~2a  , 2a*  = l-a2 

(3.5) 

a 

= EcosX,  a*  = E cos  Y 

(3.6) 

a2 

= E cos  2X,  a 2 = E cos  2Y  . 

(3.7) 

On  simplifying,  we  obtain 

Zl-2  = diag(a1ai'  aia2'  a2a2'  a2ai) 


(3.8) 
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Note  that 


A A A 

a = R^  a2  = 2 { 1 - 

- (RjA^} 

(3.9) 

* * _ ^ * _2 
a = R2,  a2  = 2 { 1 - R2 

- (R2A2n 

(3.10) 

e(U')  = (aa*,  0,  0,  0) 

(3.11) 

and 

E(T')  = (a,  0,  a*.  0) 

• 

(3.12) 

Substituting 

these  results  in  (2.3), 

we  obtain 

2 2 

2 2 

2 

r‘°i-VV  , °2 

u.  u. 

3 i 4 l 

(3.13) 

r XY 

‘ 1 Vt  VS 

a a*  ""  a*a„J  ' 
2 2 12 

where  XQ  = YQ  = 0 . 

Note  that 

a ^ 

= var(cosX)  , a2  = var(sinX) 

A 

A 

"'a*  = var(cosY)  , a*  = var(sinY)  . 

Since  (U',T')  is  jointly  normal  for  large  n , ng(X,Y) 
is  asymptotically  x^  under  H . Now  since  (2.1)  and  (2.2) 
are  asymptotically  equivalent,  it  follows  that  under  H , 


is  asymptotically 


4 . Properties  of  the  Correlation  Coefficient. 


(a)  It  can  easily  be  seen  that  the  population  counterpart 
of  (3.14)  is 


where 


(4.1) 


01  = corrjcos  (X  - uQ)  , cos(Y- vQ]  . 

P2  = corr{cos  (X -U())  , sin(Y-vQ)}  , 

O j = corr{sin(X  - Uq)  , cos  (Y  - vQ)  } , 

P4  = corr{sin(X -uQ)  , sin(Y-vQ)]  , 

and  u _ and  are  the  mean  directions  of  X and  Y 

0 0 

respectively  (Note  that  corr(X,Y)  denotes  the  ordinary 
correlation  coefficient) . 

2 

We  have  0 < P < 4 . 

2 

If  X and  Y are  independent,  we  have  D^y  = 0 . When 

2 2 2 2 2 

p =o=p=p  =1  we  have  P = 4 . In  this  case, 

1 2 3 4 xy 

X = Y=0,  with  probability  1 . For  the  perfect  rotational 
dependence  of  the  form  . Y = (±X  + \ji)mod  2tt  , it  is  found  that 


(b)  The  exact  distribution  of  U|T  under  (1.1)  does  not 
involve  nuisance  parameter  (u^.h^)  and  (p2,x2)  because 
T is  sufficient  for  (u^#h^)  , i = l,2  . Also  the  exact 


y 


2 

distribution  of  r does  not  involve  (u.,h  ) 

XY  ix 

by  the  same  argument  as  in  Cox  (1975) . 


i = 1,2 


(c)  Let  x and  Y take  values  on  a p-dimensional  torus, 
where  X'X  = Y'Y  = i . Then  (4.11)  immediately  extends  to 

P P 

pI,  = £ * {Corr(X*  Y*)}2  , 

i=l  j=l  1 1 

where  X*  = AX,  Y*=  BY  , and,  A and  B are  orthogonal 
matrices  such  that  E(X*)  = E(Y*)  = e^  , where  = ( 1 , 0 , 0,  . . . , 0) 

Its  sample  counterpart  can  be  written,  and  the  above  discussion 
for  p = 2 can  be  generalized. 


5 . Numerical  Example 

Downs  (1974)  considered  the  following  data  related  to 
the  estimated  peak  times  (converted  into  angles)  for  two 
successive  measurements  of  diastolic  blood  pressures. 


o = 30°  , 15°,  11°,  4°,  348°,  347°,  341°,  333°,  332°,  28' 
cp  = 25°,  5°,  349°,  350°,  340°,  347°,  345°,  331°,  329°,  28°. 


Then 


r2  =0.974,  r2  =0.213,  r2  =0.152,  r2  =0.933 
2 

Hence  r = 2.27  , 

c which  is  large.  Hence  there  is  a strong  evidence  for 

dependence  as  it  is  expected  on  practical  ground.  Note  that 

9 2 

nr  = 22.7  , and  the  1%  value  of  is  13.3. 
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